A high-order, analytically divergence-free discretization method for Darcy's problem

We develop and analyze a meshfree discretization method for Darcy’s problem. Our approximation scheme is based upon optimal recovery which leads to a collocation scheme using divergence-free positive definite kernels. Besides producing analytically incompressible flow fields, our method can be of arbitrary order, works in arbitrary space dimension and for arbitrary geometries. After deriving the scheme, we investigate the approximation error for smooth target functions and derive optimal approximation orders. Finally, we illustrate the method by giving numerical examples.